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ABSTRACT 

We report on the results of INTEGRAL observations of the neutron star low 
mass X-ray binary SAX J1810. 8-2609 during its latest active phase in August 
2007. The current outburst is the first one since 1998 and the derived luminosity 
is 1.1 — 2.6 X lO^^ergs"^ in the 20-100 keV energy range. This low outburst 
luminosity and the long-term time-average accretion rate of ~ 5 x lO^^^M© 
yr~^ suggest that SAX J1810. 8-2609 is a faint soft X-ray transient. During the 
flux increase, spectra are consistent with a thermal Comptonization model with 
a temperature plasma of /cTg ~ 23-30 keV and an optical depth of r ~ 1.2-1.5, 
independent from luminosity of the system. This is a typical low hard spectral 
state for which the X-ray emission is attributed to the upscattering of soft seed 
photons by a hot, optically thin electron plasma. During the decay, spectra 
have a different shape, the high energy tail being compatible with a single power 
law. This confirm similar behavior observed by BeppoSAX during the previous 
outburst, with absence of visible cutoff in the hard X-ray spectrum. 
INTEGRAL/ JEM-X instrument observed four X-ray bursts in Fall 2007. The 
first one has the highest peak flux 3.5 Crab in 3-25 keV) giving an upper 
limit to the distance of the source of about 5.7 kpc, for a LEdd ~ 3.8 x 10^^ erg 
s~^. The observed recurrence time of ~1.2 days and the ratio of the total energy 
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emitted in the persistent flux to that emitted in the bursts {a ~73) allow us to 
conclude that the burst fuel was composed by mixed hydrogen and helium with 
X>0.4. 

Subject headings: gamma rays: observations - radiation mechanisms: non-thermal 
- stars: individual: SAX J1810. 8-2609- stars: neutron - X-rays: binaries 

1. Introduction 

The transient X-ray source SAX J1810. 8-2609 was discovered on 1998 March 10 (Uber- 
tini et al. 1998) with the Wide Field Cameras (2-28 keV) on-board the BeppoSAX satellite. 
During the performed Galactic Bulge monitoring, a strong type 1 X-ray burst was detected 
identifying this compact object as a neutron star in a low-mass X-ray binary system. As- 
suming standard burst parameters and attributing the photospheric radius expansion to 
near-Eddington luminosity, the distance was estimated be ~ 5 kpc (Natalucci et al. 2000). 
The wide-band spectral data (0.1-200 keV) obtained later with the NFl/ BeppoSAX showed 
a hard X-ray spectrum described by a power law with photon spectral index F = 1.96 ± 0.04 
and a soft component which was compatible with blackbody radiation of temperature kT~0.5 
keV (Natalucci et al. 2000). 

During a ROSAT follow-up observation on 1998 March 24 an X-ray source (named RX 
J1810. 7-2609) was detected at a position consistent with the WFC error box (Greiner et 
al. 1998). Optical-to-infrared follow-up observations of the 10" radius ROSAT HRl X-ray 
error box revealed one variable object (R=19.5 ± 0.5 mag on March 13, R> 21.5 mag on 
1998 August 27) which was proposed as the optical/IR counterpart of RX J1810. 7-2609 and 
SAX J1810.8-2609 (Greiner et al. 1998). 

Using Chandra instruments, Jonker et al. (2004) detected the neutron star system 
in quiescence at an unabsorbed luminosity of ~ 10^^ergs~^ (assuming a distance of 4.9 
kpc). The quiescent spectrum is well-fitted with an absorbed power law with a photon index 
F = 3.3 ± 0.5 and the Galactic absorption value {Nu^gai ~ 3.3 x 10^^ cm~^) consistent with 
the value derived in outburst. 

Since 1998 this burster remained in a quiescent state. Only in August 2007 Swift 
observed a new phase of activity (Parson et al. 2007). The Swift/UVOT instrument detected 
a weak source in the white-band filter at the position of SAX J1810. 8-2609 and did not detect 
it in any other single filter (Scady et al. 2007). The source was observed on a daily basis with 
Swift using ~ Iksec exposure, starting August 6, 2007. In all observations the Swift/XRT 
0.3 — lOkeV spectrum was well fitted using an absorbed power law model with a hydrogen 
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column density of A^^^ ~ 5 x 10^^ cm~^ and a spectral index of F ~ 2 (Degenaar et al. 2007). 
After few months, the source went back to quiescent state, indeed on November 3rd and 5th, 
Swift/XKT did not detect it during two individual ~ 1.6ksec and ~ 1.9ksec observations. 

The X-ray spectra of LMXBs are usually fit with a complex model: at low energies a 
blackbody component that approximates the spectrum of an optically thick, geometrically 
thin accretion disk and/or the neutron star surface, and at higher energies a Comptonization 
component due to repeated inverse Compton scattering of the soft seed photons by hot 
electrons plasma with a thermal distribution of velocities. BeppoSAX and INTEGRAL 
results showed that the hard component can extend up to 200 keV without any appreciable 
break (Di Salvo et al. 2000, 2001, Fiocchi et al. 2006, Piraino et al 1999, laria et al. 2001, 
Tarana et al. 2006). In this paper, we study the spectral behavior of the X-ray transient 
burster SAX J1810. 8-2609, showing this behavior during the decay of the outburst. Finally, 
we report on four X-ray bursts observed by INTEGRAL/ JEM-X instrument in Fall 2007. 



2. OBSERVATIONS AND DATA ANALYSIS 

The INTEGRAL (Winkler et al. 2003) observations are divided into uninterrupted 
2000 s intervals, the so-called science windows (SCWs). Spectra and ligh t curves of the sou rce 



are obtained using data from the two high-e nergy instruments JE M-Xl ( iLund et al. 20031 ) in 



the 3-20 keV band and from IBIS/ISGRI (lUbertini et al. 20031 ) in the range 22 - 200 keV. 



The instrument data are extracted for each individual SCW and processed using the Off-line 
Scientific Analysis (OSA v7.0) software released by the INTEGRAL Scientific Data Centre 
(Courvoisier et al. 2003). Following the standard analysis, we use the latest response matrix 
with 64 channels. Then, data above 90 keV are rebinned to improve the signal to noise ratio. 
The RXTE/ ASM (Levine et al. 1996) daily averaged light curve, provided by the ASM/RXTE 
teams at MIT and at the RXTE SOF and GOF at NASA's GSFC, (Fig. 1, panel a, from 
http://xte.mit.edu/ASMJc.html) shows that SAX J1810. 8-2609 has been continuously ac- 
tive since beginning of August for two months with multiple peaks. The outburst of this 
transient source was frequently observed by INTEGRAL (Haymoz et al. 2007, Galis et al. 
2007) during the Key Programme on the Galactic Centre and private Target of Opportunity 
observations. 

The IBIS/ISGRI light curve (Fig. 1, panel b and c) shows a gradual brightening in two 
energy band, 22-45 keV and 45-68 keV, while the ASM peak intensity was not monitored 
with INTEGRAL . 



We report here on the outburst emission measured by the IBIS/ISGRI instrument by 
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dividing it in four separate epochs (see Table 1). These correspond to time periods during 
which the source spectra appear quite stable, with very small or absent spectral variability 
as monitored on the time scale of a few SCWs. 

We searched simultaneous IBIS and PCA (Glasser et al. 1994) data in the XTE public 
archive but unfortunately only for epoch 1 PCA standard products are available. No 
public PCA data are available for the epoch 2, 3 and 4 . For our analysis we use the PCA 
standard products OBSID 93414-01-04-01. Data are collected in standard2 modes with a 
time resolution of 16s and 129 energy channels and from PCU 2 and PCU 4. Source and 
background spectra are generated with SAEXTRCT version 4. 2d and response files with 
the tool PCARMF vlO.l. Background rates were estimated using the epoch-5 models, as 
provided by the PCA calibration team. 

3. SPECTRAL ANALYSIS 

3.1. The wide band outburst emission 

The IBIS spectra extracted for the four epochs listed in Table 1 have been fitted with 
both a simple power law and a comptt model (Titarchuk 1994), assuming a spherical geom- 
etry for the Comptonizing region. Results are reported in Table 2. The temperature of the 
Comptonizing electrons kTg and the plasma optical depth Xp were free parameters in the fit, 
while the temperature of the soft photon Wien distribution kTg was fixed at 0.6 keV. This is 



^ http://heasarc.gsfc.nasa.gov/docs/xte/xhp_archive.html 





Instrument 


Tstart 


Tstop 


Exposure 


count s/s" 






MJD 


MJD 


ksec 




1 


IBIS 


54337 


54357 


572 


5.80±0.07 


2 


IBIS 


54358 


54362 


107 


11.6±0.2 


3 


IBIS 


54367 


54370 


121 


10.9±0.2 


4 


IBIS 


54373 


54376 


128 


10.2±0.1 


5 


PCA 


54345.49 


54345.51 


2 


75.0±0.3 



Table 1: Log of SAXJ1810.8-2609 IBIS and PCA observations. 

" Rates are in the 22-200 kev energy range for IBIS and in 3-30 keV for PCA spectra. Source 
counts are background subtracted. 
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the value observed by BeppoSAX in 1998 and previously reported by Natalucci et al. (2000). 
Spectra are well described by a simple power law after the outburst peak, while a comptt 
model is required before the outburst peak. Using a thermal Comptonization model comptt 
instead of simple power law did not give significantly better fits for epoch 3 and 4, with the 
corresponding F-test chance probabilities being 6 x 10~^ and 3 x 10~^, respectively. Instead 
this model is statistically highly significant for the epoch 1 and 2, with the corresponding 
F-test chance probabilities of 4 x 10~^ and 2 x 10~^, respectively. 

For the first period we build a spectrum in a broad energy band (3-200 keV), using 
simultaneous IBIS and PCA data. The most simple model which provides a good fit to this 
spectrum is made up of a thermal comptonized component comptt in XSPEC (Titarchuk 
1994) with a spherical geometry plus a soft component consisting of a single temperature 
blackbody and an Gaussian component for iron line. 

In the fitting procedure, a multiplicative constant has been introduced to take into account 
possible cross calibration mismatches between the soft X-ray and the INTEGRAL data; 
this constant has been found to be 1.05±0.05. Results are reported in Table 2. The iron line 
centroid is 6.3lo;4 keV, crpe < O.GkeV and equivalent width 139l5oel^. Figure [2] shows four 
spectra and the residuals with respect to the corresponding best fits. Data and models are 
shown in Figure [3|, for four epochs. 

3.2. The bursts emission 

In Fig. m and [5] we show the JEM-X 3-25 keV light curves for the four bursts. 

The start time for each burst was determined when the intensity rose to 10 % above the 
persistent intensity level. The rise time is defined as the time between the start of the burst 
and the time at which the intensity reached 90 % of the peak burst intensity, as measured 
from the 2 s bin light curve in the full 3-25 keV band. The burst duration is the approximate 
time it takes to the 3-25 keV intensity (averaged over 3 consecutive bins) to decrease back 
to the average persistent level previous to the burst start. 

The spectral analysis of the bursts is based on JEM-X data in the 3-25 keV band. 
Unfortunately, time resolved spectral analysis of such short bursts, requiring relatively high 
time resolution, leads to statistically poor results due to the little aperture of the JEM-X 
instrument. Therefore, a time- averaged spectral analysis over the first 18 s including the 
peak has been performed for each burst and every burst spectrum is well fit by a simple 
blackbody model. 

The inferred blackbody temperature, fcTbb, and apparent blackbody radius at 5.5 kpc. 
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IBIS spectra 
Epochs 1, 2, 3, 4 
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keV 
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FluX20-l00keV 


xl 

x'(d.o.f.) 


1 

2 
3 
4 


r, 00+0 05 

Z.OZ_Q Qg 

-0.06 
9 /iq+0.06 




Power Law 


3.6 
7.5 
6.3 
6.0 


1.43 [30] 
1.49 [29] 
1.13[27] 
0.92 [32] 


1 

2 
3 
4 




30lf 
231^ 
6911 
87li 


comptt 

1 9+0.4 
-■-•^-0.5 

< 0.8 

< 0.8 


3.6 
7.5 
6.3 
6.0 


0.70 [29] 
0.40 [28] 
1.02[26] 
0.81[31] 



IBIS and PCA spectra 
Epoch 1 



kTbb kTe T Fluxz-iQQkeV xl 

keV keV w-w^rg cm-^ s-i X^(d-o.f.) 

1 0.44 ±0.06 22 ±3 1.7 ±0.3 8.1 1.05[78] 



Table 2: Parameter values of spectral models fitting the outburst emission in the energy 
range 22-200keV using IBIS data and in the energy range 3-200keV using PCA plus IBIS 
data. Uncertainties are given at a 90 % confidence level. 
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Fig. 1. — a) The RXTE ASM light curve daily averaged in 2-10 keV energy band of 
SAX J1810.8-2609. b) The IBIS light curve in 22-45 keV energy band. Stars indicate 
only burst times c) The IBIS light curve in 45-68 keV energy band. 

Rhh, for every burst are reported in Table [3l Burst fiuences are obtained from the bolometric 
fluxes, Fboi, extrapolated in the 0.1-100 keV energy range and integrated over the respective 
burst durations. The peak fluxes, -Fpeak, are obtained by comparing the peak count rate of 
the 2 s bin light curves with the time-averaged count rate of the spectra medeled with a 
blackbody (see Table 3). Bolometric fluxes are extrapolated between 0.1 and 100 keV using 
XSPEC software. All uncertainties in the spectral parameters are given at a 90 % confidence 
level. 

The first burst has the highest peak flux reaching a value of ^ 500 cts/s corresponding 
to ~ 3.5 Crab (1 Crab ^ 3 x lO^^erg cm^^ s^^ between 3-25 keV). From the light curve, 
it even seems to be preceded by a precursor 24 s in advance. However, since this burst was 
observed close to the limit of the JEM-X field of view, the significance of the precursor is quite 
low and therefore its reality may be doubtful. Assuming that the main peak flux corresponds 
to the Eddington limit for a helium burst, L-^dd ~ 3.8 x 10^^ erg s~^, as empirically derived 
by Kuulkers et al. (2003), we can calculate an upper limit to the distance of the source of 
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Table 3: Analysis results of the bursts 



Dataset Burst 1 Burst 2 Burst 3 Burst 4 

Date (YYYYMMDD) 20070915 20070924 20070930 20071001 

Start time (UTC) 23:20:18 19:52:50 17:10:03 21:42:10 

^^(keV) lAt^s 2.1^°| 2.61^3 2.5^^ 

i?bb,d3.s.pc (km) 18lf G.ltl A.ltl A.btl 

xVdof 17/27 40/43 30/24 31/57 

Fboi " 4.5 ±2.3 2.8 ±0.8 2.7 ±0.7 3.0 ± 1.3 
Burst parameters 



-^peak 


9.6 ± 1.9 


5.0 ±0.6 


6.8 


±0.7 


5.5 ±1 




1.1 ±0.7 


0.7 ±0.2 


0.8 


±0.2 


0.7 ±0.3 


Rise time (±1 s) 


5 


3 


4 




7 


Duration (±2 s) 


30 


30 


30 




25 




12 ±9 


14 ±5 


12 


±4 


13 ±8 


7 (10-2) 


0.9 ±0.2 


1.0 ±0.2 


0.9 


±0.2 


0.8 ±0.2 



" Unabsorbed flux (0.1-100 keV) in units of 10 ^erg cm ^ s ^. Fluence (10 ^erg cm 2). 

r(sec) = /b/i^peak- 7 = ^pers/^peak,Max; ^pers is the persistent flux in 0.1-100 keV energy 
range previous to the time of each burst, -Fpeak,Max is the highest burst peak flux, here -Ppeak of 
Burst 1. 



about 5.7 kpc. For comparison, the theoretical value of L-^dd = 2.9 x 10^^ erg s~^, assuming 
a helium atmosphere, a canonical mass of 1.4 and 10 km radius for the neutron star 
photosphere (e.g. Lewin et al. 1993), leads to a distance of 5.0 kpc, consistent with the 
distance previously derived by Natalucci et al. (2000). The three remaining bursts are all 
weaker, reach approximately the same peak flux, and have similar decay times. From the 
detection of four bursts during the total observation time of about 928 ks elapsed on the 
source by INTEGRAL, due to the non continuous coverage of the outburst, we can estimate 
an approximate recurrence time of 2.7 days in average. Nevertheless, as the fourth burst 
occurred the day after the third burst, namely At3_4 = 102730 s = 1.2 days later, this 
interval represents a more stringent constrain on the bursting rate (see below). 

The average persistent unabsorbed flux between 0.1-100 keV, Fpcrs ~ 5 x lO^''"'^ erg 
cm-2 s"^, translates to a bolometric luminosity Lpcrs ~ 1-8 x 10^^ erg s^-*^, assuming an 
approximate distance of 5.5 kpc. This corresponds to a mass accretion rate per unit area 
equal to m = Lpers(l + z) ri~^ c-^/Aacc ~ 10^ g cm-^ (where Aacc = 47ri?^g and 
77 = GMns/(-Rns c^) ~ 0.2 is the accretion efficiency for a canonical neutron star). 
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channel energy (keV) channel energy (keV) 



Fig. 2. — Four spectra of SAX J1810. 8-2609 and the residuals with respect to the corre- 
sponding best fits, COMPTT plus blackbody component for epoch 1 . simple comptt model for 
epoch 2 and simple power law component for epochs 3 and 4. For the epoch 1 spectrum is 
obtained using PCA and IBIS data (panel 1), while for others epochs spectra are obtained 
using only IBIS data (panels 2, 3 4). 

The total energy released by the first burst was Ef^^i ~ 4 x 10^^ erg which, assuming 
complete and isotropic burning, corresponds to an ignition column y = E^^i^l+z) /AnR'^gQrmc 
ranging between y Ki 1 x 10^ g cm~^ for burning hydrogen with abundance X=0.7, and 
y 2.6 X 10^ g cm~^ for X=0 (pure hehum); here Qnuc = 1-6 + 4X MeV/nucleon is the 
nuclear energy release for a given average hydrogen fraction at ignition X, and z=0.31 is the 
appropriate gravitational redshift at the surface of a 1.4 Mq neutron star (Gumming, 2003). 
From the relation Atj-cc = y(l + z)/7h a burst recurrence time of 1.5 days is expected for 
X=0.7, and Atj-cc = 3.8 days for pure helium burning. The same calculations for the fourth 
burst with an energy release of Eb^ ~ 2.5 x 10'^^ erg lead to Atrec = 0.9 days for X=0.7, and 
Atrec = 2.3 days for X=0. The observed recurrence time seems thus most consistent with 
mixed H/He burning. 

Moreover, it is also possible to calculate the burst energetics by the ratio of the total 
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channel energy (keV) 

Fig. 3. — Four different spectral states with the comptt modeL For the epoch 1 PCA and 
IBIS spectra are fitted (red color), while for others epochs only IBIS spectra are fitted (green, 
blue and orange colors). 



energy emitted in the persistent flux to that emitted in the bursts (e.g. Galloway et al., 
2004): a = {F^,J h)At,,, ^ 73, for At,,, = At^-^, F^er. = 5 x 10' 



-in —2 
erg cm s 



and 



ff, ~ 0.7 X 10~^ erg cm~^ is the fiuence of the fourth burst. Assuming that all the accreted 
fuel is burned during the bursts, the calculated a-value from the measurable quantities is 
consistent with Qnuc = {I + z)ric^ / a {10^^ erg/g)^^ 3.2 MeV/nucleon, corresponding to an 
hydrogen fraction X=0.4. Since other bursts could have been burnt during the observation 
gaps, the calculated a value is only an upper limit and conversely the calculated value of X 
is a lower limit. We can indeed conclude that the burst fuel could be composed by mixed 
hydrogen and helium with X>0.4. 



4. DISCUSSION 



The IBIS/ISGRI observations have allowed us to follow the high energy behavior of 
SAX J1810. 8-2609 during its long and bright X-ray outburst. Light curves varied simulta- 
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-20 20 40 

Time from burst start (s) 



Fig. 4. — The first type I X-ray burst detected on September 15, 2007. Time corresponds 
to 23:20:18 (UTC). The JEM-X (3-25 keV) hght curve is shown with a time bin of 2 s. 

neously in all X-ray bands monitored (2-10 keV, 22-45 keV, 45-68 keV). The X-ray spectra 
appeared always to be comparably soft, with a photon index of F = 2.3 — 2.7. The 20- 
100 keV luminosities are in the range Lx = 1-1 — 2.6 x lO^^ergs"^ (estimated for a 
source at 5.5 kpc), which is typical for the low hard state of neutron star binaries (Barret et 
al. 2000). We estimated a fluence of 5.5 x lO^^er^f cm~^ using an average bolometric flux of 
1.6 X cm~'^ in the 0.1-500 keV energy band and an outburst duration of ~40 days 

from ASM light curve. The long term time averaged accretion rate is M ~ 5 x IO^^^Mq 
yr~^, taking into account a time interval between two outbursts of ~9.6yr. The ASM, IBIS 
and BAT □ light curves show the same behavior, so the outburst duration and the time 
between two outbursts are estimated using the ASM light curve. This time-average low 
mass accretion rate, the outburst luminosity of ~ 1 — 3 x 10'^^ erg lower than typical 
values for neutron star soft X-ray transient (~ 2 x 10'^'^ erg s^^), together with low quiescent 



•^in the 15-50 energy range, from http://heasarc.gsfc.nasa.gov/docs/swift/results/transients/weak/SAXJ1810.8- 
2609 
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Fig. 5. — The three other bursts respectively detected on September 24. 30, and October 
1st, 2007. The JEM-X (3-25 keV) hght curves are shown with a time bin of 2 s. 



luminosity (~ 1 x 10'^^ erg s^^) reported by Jonker et al. 2004, strengthens the idea of these 
authors that this source belongs to the class of faint soft X-ray transient. In addition we 
note that the low average persistent bolometric luminosity is very similar as the luminosity 
of the Ultra Compact X-ray Binaries (see, e.g., Fiocchi et al. 2008 and Falanga et al., 2008). 
However, we think difficult to derive a conclusion about the ultra compact nature of the 
system because the derived hydrogen fraction in the burst fuel of SAX J1810. 8-2609 is not 
consistent with an ultra compact source, since those are thought to accrete pure Helium 
from a white dwarf (Nelemans and Jonker, 2006). 

Spectral parameters are not correlated with the observed luminosities, but instead, they 
vary according with the rise/decay phases of the outburst. During the rise of the flux, the 
SAX J1810. 8-2609 luminosity changes by a factor 2, while there are not modifications of 
the spectral shape: the electron temperature kT^ is ~ 23-30 kcV and optical depth r of the 
plasma is ~ 1.2-1.5. This hard X-ray emission could be interpreted in the standard way, as 
produced by the upscattering of soft seed photons by a hot, optically thin electron plasma. 
During the decrease of the flux, spectra show a harder spectral shape with an optical depth 
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of the plasma lower than 0.8 and very high electrons temperatures kT^ of ~ 69-87 keV. The 
spectral parameters measured during the decay phase of the 2007 outburst agree with the 
ones found using the BeppoSAX observations (Natalucci et al. 2000), showing the same 
X-ray spectral behavior: during the decay phase of the outbursts of 1998 also no high energy 
spectral steepening was observed. 

We cannot determine whether the emission is due to a thermal or non-thermal process, 
because equally good fits are obtained either with a power law with no detectable cutoff 
below ~100 keV or with a thermal Comptonization spectrum with an electron temperature 
in excess of ~ 80-90 keV. This electron distribution could also arise from Comptonization 
by hybrid (thermal and non thermal) corona (Coppi 1999), or from the Compton cloud 
located inside the neutron stars magnctosphcrc (Titarchuk ct al. 1996), or, alternatively, the 
power-law component could be produced by Comptonization of synchrotron emission in a 
relativistic jet (Bosch- Ramon et al. 2005, Fender 2004). Up to date there are few detections 
of radio emission associated with neutron star X-ray transients and sometimes outbursts of 
soft X-ray transient were associated with strong transient radio emission (Ball et al. 1995, 
Kuulkers et al. 1999, Fender & Kuulkers 2001). A comparison with hard tails detected from 
neutron star systems and some black hole binaries could be interesting suggesting that a 
similar mechanism could originate these components. 

We acknowledge the AST financial/programmatic support via contracts ASI-IR 1/008/07/0. 
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